Trigonal warping and Berry's phase Nn in ABC-stacked multilayer graphene 



o 
o 

c 

in 



Mikito Koshino 1 and Edward McCann 2 
1 Department of Physics, Tokyo Institute of Technology, 
2-12-1 Ookayama, Meguro-ku, Tokyo 152-8551, Japan 
2 Department of Physics, Lancaster University, Lancaster, LAI J^YB, UK 

The electronic band structure of ABC-stacked multilayer graphene is studied within an effective 
mass approximation. The electron and hole bands touching at zero energy support chiral quasipar- 
ticles characterized by Berry's phase N-rr for TV-layers, generalizing the low-energy band structure 
of monolayer and bilayer graphene. We investigate the trigonal-warping deformation of the energy 
bands and show that the Lifshitz transition, in which the Fermi circle breaks up into separate parts 
at low energy, reflects Berry's phase Nit. It is particularly prominent in trilayers, TV = 3, with the 
Fermi circle breaking into three parts at a relatively large energy that is related to next-nearest- 
layer coupling. For TV = 3, we study the effects of electrostatic potentials which vary in the stacking 
direction, and find that a perpendicular electric field, as well as opening an energy gap, strongly 
enhances the trigonal- warping effect. In magnetic fields, the TV = 3 Lifshitz transition is manifested 
as a coalescence of Landau levels into triply-degenerate levels. 

PACS numbers: 71.20.-b,81.05.Uw,73.63.-b,73.43.Cd. 
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I. INTRODUCTION 

Soon after the fabrication of individual graphene flakes 
a few years ago [l[ , it was realized that low-energy quasi- 
particles in graphene are chiral, with a linear dispersion 
and degree of chirality characterized by Berry's phase tt 
in monolayer graphene [H, S] and quadratic dispersion re- 
lated to Berry's phase 2ir in bilayers [J, [5( . In addition to 
their degree of chirality, bilayers are distinguished from 
monolayers by the possibility of using doping or exter- 
nal gates to induce interlayer asymmetry that opens a 
tunable gap between the conduction and valence bands 
[1, H, 0, B H, EE El j as observed in tr ansp ort fl2l. [l3l 
and spectroscopic measurements [iH [TH [ll* [l7l. Ha. Il9j] . 

Recently, there has been experimental interest in the 
transport properties of trilayer graphene [IE HJ, [2^ |. 
It is expected that two different types of stacking or- 
der, ABA and ABC (illustrated in Fig. 1), will be real- 
ized in nature and that electronic properties will depend 
strongly on the stacking type. For ABA-stacked trilayer 
graphene, the low-energy electronic band structure con- 
sists of separate monolayer-like and bilayer-like bands 
@, 0, HE [H, HE IH, that become hybridized in the 
presence of interlayer asymmetry 0, HH . By contrast , the 
low-energy bands of ABC-stacked trilayers @, 
do not resemble those of monolayers or bilayers, but ap- 
pear to be a cubic generalization of them. Thus, there is 
a cubic dispersion relation and chirality related to Berry's 
phase 3-7T @, HEHilj and, as in bilayers, the application 
of interlayer asymmetry is predicted to open an energy 
gap in the spectrum @, EE] • 

In this paper, we show that the low-energy band struc- 
ture of ABC-stacked multilayer graphene is not just a 
straightforward generalization of that of monolayers and 
bilayers. We focus on a particular aspect of the band 
structure, trigonal warping, which plays a crucial role in 
the low-energy band structure. Trigonal warping is a de- 
formation of the Fermi circle around a degeneracy point 



[3l| , at each of two inequivalent corners of the hexagonal 
Brillouin zone that are known as K points [32[ [Fig. 1(b)]. 
In bilayer graphene, trigonal warping is enhanced by the 
interlayer coupling and leads to a Lifshitz transition [33J 
when the Fermi line about each K point is broken into 
several pockets @, [H HE H HE H, M, HI . Here, we de- 
velop an effective Hamiltonian for ABC-stacked trilayer 
graphene, to show that trigonal warping in it is both 
qualitatively and quantitatively different from that in bi- 
layers. The main contribution to trigonal warping arises 
from a different type of interlayer coupling that is miss- 
ing in bilayers and we predict that it leads to a Lifshitz 
transition at a much larger energy ~ lOmeV, which is 10 
times as large as in a bilayer. Moreover, on undergoing 
the Lifshitz transition, the Fermi surface breaks into a 
different number of pockets reflecting Berry's phase 37r 
in contrast to 2tt in bilayers. Here, we also generalize our 
approach to describe trigonal warping in general ABC- 
stacked TV-layer graphene, to show that Berry's phase Nir 
manifests itself in different characteristics of the Lifshitz 
transition. 

In the next Section, we describe the effective mass 
model of ABC-stacked trilayer graphene and the result- 
ing band structure. Then, in Section IIII1 we derive an 
effective low-energy Hamiltonian and we use it to com- 
pare the behavior of low-energy chiral quasiparticles in 
trilayers with those in monolayer and bilayer graphene. 
In Section [IV! we provide an approximate analytical de- 
scription of the Lifshitz transition in the absence and in 
the presence of interlayer asymmetry that opens a gap 
in the spectrum. Section [V] describes the manifestation 
of the Lifshitz transition in the degeneracy of Landau 
levels in the presence of a finite magnetic field. In Sec- 
tion |Vll we generalize our approach to ABC-stacked Al- 
layer graphene. Throughout, we compare the approxi- 
mate description of the effective low-energy Hamiltonian 
with numerical diagonalization of the full effective mass 
model. 
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FIG. 1: (a) Schematic of the ABC-stacked trilayer lattice 
containing six sites in the unit cell, A (white circles) and B 
(black circles) on each layer, showing the Slonczewski-Weiss- 
McClure parameterization 391 of relevant couplings 70 to 74. 

(b) Schematic of the hexagonal Brillouin zone with two in- 
equivalent valleys K± showing the momentum p measured 
from the center of valley K+. Schematic of the unit cell of 

(c) ABC-stacked trilayer graphene, (d) ABA-stacked trilayer 
graphene, and (e) bilayer graphene. In (c), 72 describes a 
vertical coupling between sites B3 and Al in different unit 
cells. 



II. THE EFFECTIVE MASS MODEL OF 
ABC-STACKED TRILAYER GRAPHENE 

The lattice of ABC-stacked trilayer graphene consists 
of three coupled layers, each with carbon atoms arranged 
on a honeycomb lattice, including pairs of incquivalcnt 
sites {A1,B1}, {A2,B2}, and {A3,B3} in the bottom, 
center, and top layers, respectively. The layers are ar- 
ranged as shown in Fig. nja,c), such that pairs of sites 
Bl and A2, and B2 and A3, lie directly above or below 
each other [for comparison, the unit cell of ABA-stacked 
graphene is shown in Fig. [TJd)]. In order to write down 
an effective mass Hamiltonian, we adapt the Slonczcwski- 
Weiss-McClure parameterization of tight-binding cou- 
plings of bulk graphite [39|]. Nearest-neighbor (Ai-Bi for 
i = {1,2,3}) coupling within each layer is described by 
parameter 70, 71 describes strong nearest-layer coupling 
between sites {B1-A2 and B2-A3) that lie directly above 
or below each other, 73 (74) describes weaker nearest- 
layer coupling between sites A1-B2 and A2-B3 {A1-A2, 
B1-B2, A2-A3, and B2-B3). With only these couplings, 
there would be a degeneracy point at each of two in- 
equivalent corners, K±, of the hexagonal Brillouin zone 
|32j | but this degeneracy is broken by next-nearest-layer 
coupling 72, between sites Al and B3 that lie on the 
same vertical line [13, [H, ■ For typical values of bulk 



ABA graphite we quote [39] 70 = 3.16eV, 71 = 0.39eV, 
72 = -0.020eV, 73 = 0.315eV and 74 = 0.044cV. Al- 
though the atomic structures of ABA and ABC (rhom- 
bohcdral) graphite are different, we refer to those values 
in the following numerical calculations, assuming that 
the corresponding coupling parameters have similar val- 
ues m. 

In a basis with atomic components ipAi, ^Pbi, 4'A2, 
2, i/Ma, *0B 3, the ABC-stacked trilayer Hamiltonian 



a abc = I d 2 v 
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where the 2x2 blocks are 
D, 



V 



Ui VTV 

vir Ui 
7i 



(*= 1,2,3), 



v 3 n 
-V4ir'< 



W 



72/2 





(1) 



(2) 



(3) 



where v = (\/3/2)a7o/?i, W3 = {^/3/2)a^z/K v a — 
(V3/2)aj4,/H, -K = £p x + ip y , = £p x - ip y , and £ = ±1 
is the valley index. Here p = (p x ,Py) — p{cos4>, sin0) is 
the momentum measured with respect to the center of 
the valley [Fig. Epa)]. The parameters Ui, U2, and U3 
describe on-site energies of the atoms on the three layers 
that may be different owing the presence of substrates, 
doping, or external gates. In the following, we set the 
average on-site energy to zero U\ + U2 + U3 = and 
write differences between the on-site energies in terms of 
asymmetry parameters Ai and A2 |26j |. 

Ai = {U x -U z )/2, 

A 2 = {U x -2U 2 + U 3 )/&. 

Parameter Ai describes a possible asymmetry between 
the energies of the outer layers, whereas A2 takes into 
account the possibility that the energy of the central layer 
may differ from the average outer layer energy. 

As there are six atoms in the unit cell, ABC-stacked 
trilayer graphene has six electronic bands at low en- 
ergy as plotted in Fig. 2. For no intcrlayer asymme- 
try Ai = A2 = 0, and exactly at the K point, p = 0, 
the eigenvalues e of Hamiltonian Eq. |T]) are given by 
(e 2 — 7 2 ) 2 (e 2 — 7|/4) = 0. Four of the bands are split 
away from zero energy by interlayer coupling 71 (e = ±71 
twice). These high-energy bands correspond to dimcr 
states formed primarily from orbitals on the atomic sites 
Bl, A2, B2, and A3 that are strongly coupled by 71. 
The other two bands (e = ±72/2) are split slightly away 
from zero energy by next-nearest layer coupling 72/2 that 
connects atomic sites Al and B3 [nil . I23I l28j . 

Figure [2] shows the band structure at several Ai's with 
A 2 = 0, using the parameter values quoted above. Ai 
opens an energy gap between the lower electron and hole 
bands, because of the energy difference between Al and 





vp x / y. 



FIG. 2: Band dispersion of ABC-stacked trilayer graphene 
in the vicinity of K+ along p x axis. Parameter values are 
70 = 3.16eV, 7i = 0.39eV, 72 = -0.020eV, 73 = 0.315eV and 
74 = 0.044eV B. 



B3 sites 0)[l3|- Figure [3] shows contour plots of the lower 
electron band at (a) A1/71 = and (b) 0.4, showing that 
the band is trigonally warped, and the contour splits into 
three pockets at low energy. The detailed band structure 
and its relation to the band parameters will be studied 
in the following sections. 
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III. THE LOW-ENERGY EFFECTIVE 
HAMILTONIAN 



FIG. 3: (a) Equi-energy contour plots of the lowest electron 
band of ABC trilayer graphene at (a) Ai = and (b) O.471. 
Numbers on the contours indicate energy in units of 71. Filled 
and empty triangles represent local minima and maxima, re- 
spectively, of the energy band. 



To describe the low-energy electronic properties of 
ABC-stacked trilayer graphene it is useful to derive an 
effective two-component Hamiltonian that describes hop- 
ping between atomic sites AI and S3. Such a procedure 
has been applied to bilayer graphene [5[ and to ABC- 
stacked trilayer graphene [7] for 73 = 72 = A 2 = 0. We 
begin with the energy eigenvalue equation H^> — e^P of 
the six-component Hamiltonian, Eq. ([T]), eliminate the 

dimer components \ — {' l pBi,ipA2,ipB2,'>pA3) T and, then, 
simplify the expressions for low-energy components 9 = 
(ipAi, r 4 , Bz) T by treating interlayer coupling 71 as a large 
energy scale \e\,vp, \y 2 \, ha\, |Ai|, |A 2 | < 71. We de- 
note hg as the diagonal block of Hamiltonian of Eq. ([1]) 
corresponding to 6, h x as the four by four diagonal block 
corresponding to x, and u as the off-diagonal 2x4 block 
coupling 9 and x- The Schrodinger equation for 9 can be 
expanded up to first order in e as [h e - u/i~V]0 = eS9 
with S = l + w/i~ 2 ut . Then, the effective Hamiltonian for 
9 = S 1 ' 2 9 becomes » S-^ 2 [h e 

Thus, we find the following two-component Hamilto- 



uh^u^S- 1 ' 2 . 



nian in a basis of the AI-B3 sites: 



H 



fflc = H 3 + H 3w +H 3c + H A1 +H A2 , (4) 
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Here we keep only the leading order for the terms includ- 
ing 79, v 3 and V4. Terms H 3 and H A i were derived in 
Ref. [3]. The cubic term H 3 describes effective hopping 
between sites ^41 and B3 via the other sites on the lattice 
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that are strongly coupled by 71. Taken on its own, it pro- 
duces a dispersion e = ±v 3 p 3 /j 2 . H 3w arises from the 
skewed interlayer coupling 73 and the next-nearest inter- 
layer coupling 72, and is responsible for trigonal warping 
as discussed in detail later. H^,,, coming from another 
interlayer coupling 74, gives an identical curvature to the 
electron and hole bands and thus introduces electron- 
hole asymmetry. Terms Hai, Ha2 arise from the inter- 
layer asymmetries Ai and A 2 , respectively. Hai leads 
to the opening of an energy gap between the conduc- 
tion and valence bands, while Ha 2 produces electron- 
hole asymmetry in a similar way as H^ c . In the two- 
component basis of H^bc, time reversal is described by 
H* (p, Ax,£) = H (— p, Ai, — £) and spatial inversion by 
<r x H(p, A!,^)a x = H (-p, -Ai,-f). Manes et al [H 
showed that the Fermi points of ABC-stacked multilayers 
are stable with respect to the opening of a gap against 
perturbations that respect combined time reversal and 
spatial inversion, as well as translation invariance. 

The low-energy effective Hamiltonian for ABC-stacked 
trilayer graphene bears some resemblance to that of bi- 
layer graphene Q. In the lattice of bilayer graphene, 
Fig. 1(e), two of the sites (Bl and A2) are directly above 
or below each other and are strongly coupled by inter- 
layer coupling 71 whereas two sites (Al and B2) do not 
have a counterpart in the other layer. The low-energy 
Hamiltonian is written in a basis (tpAi, ipB2) of these two 
sites: 

h ab = H 2 +H 2w + H A , (5) 

H 2w = v 3 ( ° t I ) , 
~ _ 2vv 4 p 2 ( 1 \ 

H2c ~ 1)' 

*>- A H£)(l-°0- 

where parameter A describes interlayer asymmetry be- 
tween on-site energy A of the atoms, Al and Bl, on 
the first layer and —A of the atoms, A2 and B2, on the 
second layer. 

The first term in each Hamiltonian, H 2 for bilayers, 
Eq. (O, and H 3 for ABC-stacked trilayers, Eq. ([4]), are 
members of a family of Hamiltonians Hj = F(p) cr ■ n 
where n = l x cos(J£(j))+ly sin( J£</>) for p = p(cos </>, sin 4>) 
[E 0, HI, H3|- They describe chiral quasiparticles, and 
the degree of chirality is J — 1 in monolayer graphene, 
J = 2 in a bilayer, and, here, J = 3 in ABC-stacked 
trilayer. Quasiparticles described by the Hamiltonians 
Hj acquire a Berry's phase — i § c dp • (\I r |Vp|\I f ) = J£,n, 
upon an adiabatic propagation along an equi-energetic 
line C. Thus charge carriers in ABC-stacked trilayer 
graphene are Berry's phase 3£7r quasiparticles, in con- 
trast to Berry's phase £-7r particles in monolayers, 2£tt 



in bilayers. As well as the first term in the Hamilto- 
nian Eq. of ABC-trilayers being a generalization of 
that in bilayers, the influence of interlayer asymmetry 
Ai = (Ui — Us)/2 as described by Hai is similar to that 
in bilayers as described by Ha, Eq. ([5]). 

IV. TRIGONAL WARPING AND THE 
LIFSHITZ TRANSITION 

In a similar way to bulk graphite 39] , the parameter 73 
[where v$ = (^/i/2)a"f 3 /K\ produces trigonal warping in 
bilayer graphene [5| , where the equi-energetic line around 
each valley is stretched in three directions. This is due 
to the interference of the matrix elements connecting Al 
and B2, where an electron hopping from Al to B2 ac- 
quires a factor e 2% ^ in H 2 and e~ 1 ^ in H 2w . We neglect 
the terms including v 4 which add a term ex p 2 to the en- 
ergy but don't contribute to trigonal warping. At A = 0, 
the eigenenergy of Eq. (J5J) is given by 

~ J- / 2 2 OC V 3V 2 P 3 Q , . V 4 p 4 

ewiWwyr — 24 cos 30 H p. (6) 

The warping has a dramatic effect when H 2 and H 2w 
have comparable amplitudes, i.e., v 2 p 2 ~ v^p, which 
is satisfied at y ~ p n = _^\ v ^ly 2 - ft leads to a Lifshitz 
transition 0, [23, HE EE [ll, [S3, [H, [13, HI] , in which the 
equi-energetic line is broken into four separate pockets. 
There is one central pocket located around p = and, 
three "leg" pockets centered at momentum of magnitude 
p = po and angle 4>o = 2mr/3 + (1 — The Fermi 

pocket separation occurs at energy ~ (v 3 /v) 2r )i/ A, 
which is estimated to be ~ lmcV. 

In ABC-stacked trilayer graphene, there is a similar, 
but much greater warping effect. In hopping from Al to 
B3, an electron acquires a factor e 31 ^ from H3 and a 
factor of unity from H$ w , giving trigonal symmetry in <j). 
At Ai — A 2 = 0, the eigenenergy of Eq. ([4} reads, 

e « ±y/ f(p) 2 + 2£/(p) g(p) cos 30 + g(pf , (7) 

where f(p) = v 3 p 3 /jf comes from H3 and g(p) — 
—2vv3p 2 /ji + 72/2 from H 3w . Similarly to the bilayer, 
the warping effect is prominent when \g(p)\ ~ f(p), or 
p ~ po with wpo/71 = [72/(27i)] 1/3 ~ («3/«)/3. This es- 
timate is valid as long as |f3/f | ^l^/Ti 1 , which holds 
for typical parameter values of bulk graphite [3{| . 

The major difference from the bilayer is the contribu- 
tion of the parameter 72/2, which appears in the Hamil- 
tonian without an accompanying momentum-dependent 
factor and, thus, it doesn't vanish at p = 0. Such trigonal 
warping produces a Lifshitz transition at low energy, but, 
unlike bilayers, it occurs at energy w I72/2I. Although 
the value of 72 in ABC-trilayer graphene has not been 
measured experimentally, comparison with similar cou- 
plings in bulk graphite [39J suggest that I72I ~ 20meV. 
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This opens up the possibility that the Lifshitz transi- 
tion in ABC-trilayer graphene could occur at a much 
higher energy than that in bilayers. At energy lower than 
I72I/2, the contour splits into three leg pockets centered 
at p ~ po in a trigonal manner. Unlike bilayer graphene, 
the central pocket is missing because H l2 does not vanish 
at p = 0. 

An effective Hamiltonian in the vicinity of the leg pock- 
ets, for |e| <§; el, may be obtained by transforming to 
momentum q = (q x ,q y ) measured from their centers, 

q x = p x cos 0o + Py sin O - Po , (8) 
q y = ~Px sin (/>o + p y cos O , (9) 

and taking the limit of infinitely large 71: 



ABC — JU 



72 


2/3^ 


271 







£aq x + iq y 



Haq x -iq v ) ( . |()) 



where a = 1 + {Av 3 /3v){2-f 1 /j 2 ) 1 / 3 . Thus, 
the pockets are elliptical with dispersion e sa 
±3|7 2 /(27i)| 2 / 3 wy / a 2 q2 + q 2. The different nature of the 

Lifshitz transition in bilayers and ABC-stacked trilayers 
is a manifestation of Berry's phase. In trilayer graphene, 
the geometrical phase integrated around the equi-energy 
line of each pocket is £7r as in a monolayer, giving 3^7r 
in total. This is different from bilayers, where 3£7r arises 
from three leg pockets and — £7r from the center pocket 
gives 2£tt in total [HUlt 

Interlayer asymmetry Ai opens a gap in the spectrum 
and produces a Mexican hat feature in the low-energy 
dispersion. Q The eigenenergy corresponding to Eq. (|4]) 
is given by 

e « ±Jf(p) 2 + 2U(p) g( P ) cos 30 + g{pf + h{p)\ (11) 



with an extra term as compared to Eq. ([7]), h(p) = 
Ai (1— f 2 p 2 /7 2 ), coming from Ha 1 ■ For no trigonal warp- 
ing [g(p) = 0], it yields e 2 = A 2 (l - v 2 p 2 /j 2 ) 2 + v 6 p 6 /jf. 
The energy is e = ±Ai at zero momentum, but there is a 
minima located isotropically about the center of the val- 
ley at finite momentum p — pi w (2/3) 1 / 4 -\/| A171 \/v 
(for I Ai| <C |7i|) at which the energy is e = e\ rs 
±A 1 (1-(2/3) 3 / 2 |A 1 / 7 i|). 

In the presence of trigonal warping, there is an in- 
terplay between the Mexican hat feature and the Lif- 
shitz transition. In the large gap regime, such that 
\g\ <C f,h, the circular edge of the band bottom is 
trigonally distorted by the perturbation of <?(p), mak- 
ing three pockets on it. The bottom of the pockets 
moves to momentum p — Pi + Spi with vSpi/ji ~ 
(v / 6/8)[7 2 /(2A 1 )] - (5/6)(« 3 /w), and energy e = e x - 
Sex with ei = (2/3) 3 /V A i/7i 



72 /2 - v ^73(« 3 /w)A 1 

The area of the pocket in k-space becomes of the order 
of PiSpi, and the depth in energy is of order Se\, both of 
which increase as increases. This significant enlarge- 
ment of the trigonal pockets, in the presence of finite Ai, 



is illustrated in Figure [3] which is produced by numerical 
diagonalization of the full Hamiltonian Eq. |T]). Note that 
similar widening of the pockets by the gap term occurs in 
bilayer graphene as well. This can be understood in an 



analogous way, by writing f(p) 
and h{jp) = A(l - 2v 2 p 2 /jf). 



v P 111, g{p) = v 3 p, 



V. LANDAU LEVEL SPECTRUM 

The energy levels in a magnetic field are given by 
replacing p with p + eA in the Hamiltonian Eq. |T]), 
where A(r) is the vector potential corresponding to the 
magnetic field. Here we consider a uniform magnetic 
field B applied along +z direction in a Landau gauge 
A = (0, Bx). Operators 7r and tt^ are then related 
to raising and lowering operators and a of the Lan- 
dau level in a conventional two-dimensional system, such 
that [l B /(V2h)]^ = at and a for K+ and K- 7 respec- 
tively, with Ib = \Jh/{eB). The operator a acts as 
atp n ,k = */nipn-i,k, and atp = 0, where tp n , k (x,y) oc 
e lky e~ z t 2 H n {z) is the wavefunction of the nth Lan- 
dau level in a conventional two-dimensional system with 
z = (x + kl 2 B )/lB, and H n being a Hermite polynomial. 

In the simplest model including only 70 and 71 without 
trigonal warping, the effective Hamiltonian H3 in Eq. (j4]) 
yields the eigenstates for K+ [3] 



e n = 0, ^ nk oc 



<Pn,k 





e«n = s—2-yJn(n - l)(n - 2) 

<Pn,k 
Stp n -3,k 



^snk OC 



(n = 0,l,2), 

> (n>3), (12) 



where s = ±1 desc ribes the electron and hole levels, re- 
spectively, As = v2hv 2 eB. The eigenstates n = 0, 1,2 
have a non-zero amplitude only on the first element {AI), 
and remain at zero energy regardless of the magnetic field 
strength, while the energy of the other levels behaves as 
oc £? 3 / 2 . At the other valley K_ , there is a similar struc- 
ture except that the first and second elements are inter- 
changed, i.e., the zero-energy Landau levels have ampli- 
tudes only on sites B3 @- 

Trigonal warping gives a remarkable feature in the 
structure of Landau levels. In enough small fields, 
the three leg pockets independently accommodate an 
equal number of Landau levels so that they are triply 
degenerate. This is in contrast to bilayer graphene 
where the central pocket also contributes to the de- 
generacy The low-energy effective Hamiltonian, 
Eq. (fT0|) . shows that the Landau level energy follows a 
similar sequence as that in monolayer graphene, e„ — 
3|7 2 /(27i)| 2/ ' 3 y / (iABSgn(n)v / n where n is integer. The 
total number of Landau levels accommodated in each 
pocket is roughly estimated by the condition e„ ~ I72I/2, 
asn~( 7l /A B ) 2 [7 2 /(2 7l )] 2 / 3 /(9a). 
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Fig- Eta) shows the Landau level spectrum at the valley 
K+ as a function of As(oc \/~B), numerically calculated 
for the full parameter model Eq. (JXJ) at Ai = A2 = 0. 
Below e = 72/2, the Landau levels are triply degenerate 
and move in proportion to \[B. The degeneracy of each 
level is broken at e — 72/2, and it splits into three sepa- 
rate levels, corresponding to coalescence of the leg pock- 
ets at the Lifshitz transition. At even higher energy, it 
approaches B 3 / 2 behavior as described in Eq. (fT2"|) . The 
triply degenerate level around zero energy is regarded 
as the n = level in each of three pockets. In actual 
fact, its degeneracy is split slightly in a large magnetic 
field, owing to magnetic break down among the semi- 
classical orbits in the leg pockets, which is caused by the 
parameter W4. When the trigonal warping vanishes, those 
three levels switch to the degenerate levels with indices 
n = 0,1,2 in Eq. (fT2"|) . 

Fig. 0Jb) shows the Landau level spectrum at K + as 
a function of asymmetry Ai with fixed magnetic field 
Ab = O.I71 (B ~ IT). As Ai is changed from negative 
to positive, three Landau levels [indicated by the single 
diagonal line that crosses e = at Ai = in Fig. EJb)] 
are pumped from the hole side to the electron side. In the 
approximate model of Eq. (fT2"|) , this corresponds to the 
fact that the energy levels n = 0, 1, 2 have a wave ampli- 
tude only on AI, so that it acquires on-site energy +Ai 
in the first order of perturbation. At the other valley 
there is the opposite movement, i.e., the three levels 
go down from positive to negative energies in increasing 
Ax. 

The energy of the Lifshitz transition appears as a re- 
gion where the levels are densely populated, and below 
that energy the levels are triply degenerated [indicated by 
the shaded region in Fig. QIb)] . It should be noted that 
the number of triply-degenerate levels increases for larger 
Ai, reflecting the enlargement of the trigonal pockets 
discussed above. In a measurement of Hall conductivity, 
those triply-degenerate Landau levels would be observed 
as quantum Hall steps of magnitude 3g v g s e 2 /h, where 
9v = 9s = 2 are the valley and spin degeneracies, respec- 
tively. 
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VI. GENERAL ABC-STACKED MULTILAYER 
GRAPHENE 



The analysis of ABC-stacked trilayer graphene can be 
extended to multilayers with N layers. We consider each 
layer to consist of carbon atoms on a honeycomb lattice, 
and the layers are arranged with ABC stacking. The 
Hamiltonian is written in a basis ipAi, "4>bi, "4>A2, ^-82, 
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FIG. 4: Landau levels of ABC trilayer graphene, plotted 
against (a) B 1 ^ 2 at fixed Ai = 0, and (b) Ai at fixed mag- 
netic field (2hv 2 eB) 1/2 = O.I71 (B ~ IT). The region in which 
Landau levels are triply degenerate is highlighted by shading. 



ipAN, 4>bn, as [40|,|41j 



H 



N 



(Di V W 

D 2 V W 

V 



(13) 



where the 2 x 2 blocks Di, V, W are defined in Eqs. (|2"l3"l) . 
Pairs of sites B{i) and A(i + 1) (i = 1, ■ ■ ■ , N — 1) are 
vertically above or below each other, and are strongly 
coupled by 71 giving dimer states. Thus, all the sites 
in the lattice, except two, contribute to bands that lie 
away from zero energy. The remaining two sites, A\ 
and Bn , form the lowest-energy electron and hole bands. 
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Note that these sites lie on the outer layers, so that the 
lowest bands are missing in an infinite system with pe- 
riodic boundary conditions applied in the stacking di- 
rection. The band structure has trigonal symmetry for 
any N. This is checked by applying the transformation 
<f> — ► cf>+2ir/3 to Eq. (Tl3|) . where the change in the matrix 
elements can be canceled by the gauge transformation 
4>An = a n ipAn and ^ Bn = a n ip Bn , with a n = e ^ 2n7r / 3 . 

The effective low-energy Hamiltonian is obtained by 
treating terms other than 7! as perturbations. The ef- 
fective Hamiltonian in a basis {4>ai, ^bn} reads 



0.04 



N ~ 

X(p) = 



\ 2vv 4 p 2 




9 


1 7l 2 


G 





X(p) 

(ni+n 2 — n ' 

{n 1 ,n 2 ,n 3 } 



1 



«3 



(vpe^r(v 3 pe-^r (y)" 3 , (14) 



where the summation is taken over positive integers 
which satisfy n\ + 2ni + 3n 3 = N. Here we collected 
all the higher order terms not including U4, but retain 
just the leading term for V4. The trigonal warping struc- 
ture can be described well in this treatment as shown 
below, since v 4 only gives the circularly-symmetric band 
curvature as in ABC trilayer. 

The eigenenergies are given by e = 2vv4p 2 /jf ± 
|X(p)|. If we neglect 72 and V3, we have X — 
(vpe^) N /(— 7i) A '~ 1 which gives a pair of bands, 
isotropic in momentum, which touch at the origin 0, 0, 
l29l . Ho]. Berry's phase integrated along an energy con- 
tour is N^n at every energy. Perturbation by 72 and 
i>3 produces trigonal warping as observed in the trilayer. 
Figure [5] shows the lower energy band structure for = 
at several iV's, where the solid lines are calculated using 
the original Hamiltonian Eq. (fT3| . and the dashed lines 
use Eq. (|14[) . We can see that the effective Hamilto- 
nian reproduces the original band structure rather well 
including the positions of the band touching points, ex- 
cept that the magnitude in energy tends to be overesti- 
mated around vp ~ 71 where the perturbative approach 
fails. 

The band touching points, or Dirac points, are given by 
the solution of X(p) = 0. They appear in a series of p's 
at only three angles fa — 2rt7r/3+(l — £)n/6, and around 
which the Hamiltonian has a chiral structure similar to 
monolayer graphene. We empirically found that the ar- 
rangement of these points obeys the following rules: We 
have [(N + l)/3] Dirac points at p 7^ at each of three 
angles, and each of them has Berry's phase £ir. Here 
[x] represents the greatest integer which does not exceed 
x. The Dirac point at the center (p — 0) only appears 
when N is not a multiple of 3, and its Berry's phase is 
£-7r and — £-7r when N = 1 and —1 (mod 3) respectively. 
The total Berry's phase summed over all Dirac points 
is always N£tt, the same as the value without trigonal 
warping. The energy scale for fine structure around the 
Dirac points becomes smaller as N increases, because the 



-0.02 



-0.04 



0.04 



0.02 



-0.02 



-0.04 



0.04 



0.02 



-0.02 



-0.04 



' / 




1 1 1 < 1 
N=2 


w 


F 


> 

y 






1 " P» 









0.5 





0.5 



FIG. 5: Low-energy band structure of ABC-stacked mul- 
tilayer graphene for several different layer numbers N, at 
Ai = A2 = 0. Solid and dashed curves are calculated us- 
ing Eq. (|13p and its approximation Eq. (|14p . respectively. In- 
sets show the equi-energy lines at e = O.O4.71. The black and 
white arrows (circles in insets) represent Dirac points having 
Berry's phase £n and — £,ir, respectively. 



matrix elements connecting A\ and Bn become higher 
order in p for larger N. We see that TV = 3 has the most 
prominent structure, where 72 directly connects A\ and 
Bn. The parameter W4 never opens a gap at the Dirac 
points but gives an energy shift by 2vv4p 2 /jf and as- 
sociated band curvature, leading to misalignment of the 
Dirac point energies as shown in Fig.[5j The curvature is 
independent of N because it is due to the second order 
process from A\ or Bn to the nearest-neighboring dimcr 
state. 

The approach applied to the Landau levels of the tri- 
layer in Sec. [V] can be extended to the A^-layer case. In 
the simplest model including only 70 and 71 , the Landau 
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levels at K + read 



e n = 0, oc 



fn,k 



(n = 0,l,---,7V-l), (15) 



A* 



Csn ~~ ° (N-l) 
^ snk I 



jry/nin - 1) • • ■ (n - N + 1) 



> (n>N), 
(16) 



with s = ±1. The first and second elements are again 
interchanged at the other valley K~. The zero-energy 
level is now TV-fold degenerate per valley and per spin 
0, US H3] • In presence of trigonal warping and V4 , how- 
ever, this is expected to split in accordance with the dis- 
crepancy between the energies of different Dirac points 
shown in Fig. [5] for B = 0, while some levels keep three- 
fold degeneracy owing to trigonal symmetry as in the tri- 
layer case. It is possible that electronic interactions may 
create exotic collective modes in such highly-degenerate 
Landau levels, but we leave the discussion of this for fu- 
ture studies. 



VII. CONCLUSIONS 

In ABC-stacked multilayer graphene with N layers, 
two low-energy bands in the vicinity of each valley are 



formed from two electronic orbitals that lie on the bot- 
tom and top layers of the system. Such bands support 
chiral q uasiparticles corresponding to Berry's phase Ntt 
[E S HO] • The interplay between different types of 
interlayer coupling produces trigonal warping, in which 
the Fermi circle around each valley is stretched in three 
directions. At very low energy, trigonal warping leads 
to a Lifshitz transition [33| when the Fermi circle breaks 
up into separate pockets, in such a way that the total 
Berry's phase is conserved. We predict that the Lifshitz 
transition is particularly prominent in trilayers, N = 3, 
with the Fermi circle breaking into three parts at a rel- 
atively large energy that is related to next-nearest-layer 
coupling. 
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